Journal of Superconductivity, Vol. XXX, 1998 



Tunneling and Photoemission in an SO (6) Supercon- 
ductor 



R.S. Markiewicz" C. Kusko" " and M.T. Vaughn" 



Received 1 June 1998 



Combining the results of tunneling, photoemission and thermodynamic studies, the pseudogap 
is unambiguously demonstrated to be caused by Van Hove nesting: a splitting of the density of 
states peak at (tt, 0). The fact that the splitting remains symmetric about the Fermi level over 
an extended doping range indicates that the Van Hove singularity is pinned to the Fermi level. 
Despite these positive results, an ambiguity remains as to what instability causes the pseudo- 
gap. Charge or spin density waves, superconducting fluctuations, and flux phases all remain 
viable possibilities. This ambiguity arises because the instabilities of the two-dimensional Van 
Hove singularity are associated with an approximate SO(6) symmetry group, which contains 
Zhang's S0(5) as a subgroup. It has two 6-component superspins, one of which mixes Zhang's 
(spin-density wave plus d-wave superconductivity) superspin with a flux phase operator. 
This is the smallest group which can explain striped phases in the cuprates. 
Evidence for a prefered hole density in the charged stripes is discussed. 
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1. 11 YEARS of the VAN HOVE SCE- 
NARIO: REMINISCENCES of R.S.M. 

Early work on the Van Hove model of high-Tc 
superconductivity was criticized for a number of rea- 
sons, chiefly (1) the perfect nesting of the Fermi sur- 
face would cause instabilities which would overwhelm 
superconductivity; and (2) the Van Hove singularity 
(VHS) is associated with a special doping, and hence 
would require 'fine tuning' of the model parameters 
to fall at the Fermi level, whereas superconductiv- 
ity exists over an extended doping range. In June, 
1987, I suggested[l] the picture which has become 
the heart of the generalized Van Hove scenario [2]. 
First, introduction of next-nearest-neighbor hopping 
~ i' (too) in a one (three)-band model - would push 
the VHS off of half filling, leading to imperfect nest- 
ing, with residual hole pockets and ghost Fermi sur- 
faces [3]. Secondly, this very sensitivity to instability 
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could in turn lower the free energy precisely when the 
Fermi surface coincides with the VHS, thereby sta- 
bilizing this special doping - and that the reason su- 
perconductivity might persist over such an extended 
doping range is that the material is inhomogeneous, 
with one phase pinned at the VHS. Due to charging 
effects, such phase separation would be nanoscale[4]. 

While a microscopic model only appeared in 
1989 [5], I quickly found that the doping dependence 
of both resistivity and Hall effect could be under- 
stood in the context of a percolation model[3,6], and 
Bill Giessen and I noted that the Uemura plot [7] 
could be interpreted in terms of an optimal doping 
for each cuprate, which we identified with the pure 
VHS phasejs]. At the 1988 MRS Meeting, Jim Jor- 
gensen asked me if the phase separation really had 
to be nanoscale, and I said it was a problem of the 
counterions: if the Sr in La2-a;Sra;Cu04 could be 
made mobile, there was no reason that a macroscopic 
phase separation couldn't occur - with one phase 
pinned at optimal doping. I didn't know then that 
Jim and co-workers had just discovered such a phase 
in La2Cu04+5, where the doping is due to mobile 
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interstitial oxygens [9]. 

One problem was: just where is the VHS. A VHS 
at half filling could produce an antiferromagnetic in- 
stability with large pseudogap[10] but no hole pock- 
ets [11], while a VHS at finite doping could explain 
hole pockets and optimal doping, but not antiferro- 
magnetism. I explored this issue by putting t' into 
a slave boson model, to separate Mott and nesting 
instabilities [5]. I found that the Mott instability - 
and hence the accompanying antiferromagnetism - 
remain locked to half filling, while for the VHS I 
found a striking surprise: correlation effects pin the 
Fermi level close to the VHS over an extended dop- 
ing range, extending across the undcrdopcd regime 
from half filling to optimal doping (close to the bare 
VHS). Numerous subsequent studies have confirmed 
the Van Hove pinning (cited on p. 1223 of Ref. [2]). 

In this same paper [5], I demonstrated a viable 
mechanism for nanoscale phase separation: a charge- 
density wave (CDW) instability is strongly enhanced 
when the Fermi level coincides with the bare VHS, 
leading to a free energy minimum at optimal dop- 
ing. Since the Mott instability produces a separate 
minimum at half filling, the energy is lowered by sep- 
arating into the two end phases. This model was 
clearly distinguished from other models, in that there 
is a highly unusual and characteristic fractional hole 
doping of the hole rich phase (this doping fixed by 
the VHS), whereas all other models had a doped 
phase with one hole per Cu. Now I had an embarras 
de richesses: both pinning of the VHS by correla- 
tion effects and phase separation, which amounted 
to a different pinning mechanism. Put another way: 
were the cuprates characterized by phase separation 
or by the (pseudo)gap associated with nesting? I 
puzzled over this issue for a number of years, con- 
vinced that somehow both viewpoints were correct 
- rather like the particle vs wave problem in quan- 
tum mechanics. I was (and am) convinced that 
both the low-temperature orthorhombic (LTO) and 
low-temperature tetragonal (LTT) phases in LSCO 
was driven by this mechanism, and I introduced 
the concept of Van Hove-Jahn- Teller effect to ex- 
plain the LTO phase as a dynamic Jahn- Teller (JT) 
phase [12, 13]. But there was something missing: there 
should have been a tetragonal dynamic JT phase at 
higher temperatures. In the years from 1987-1991, I 
was often asked where was the evidence for a CDW 
or other nesting instability. 

Things changed when I realized that the 'spin 
gap' and thermodynamic[14] data were consistent 
with a simple model for the pinned Van Hove phase 



[12,15]. New photoemission data[16] showed that 

the pseudogap is associated with the band dispersion 
near (7r,0) - i.e., the locus of the VHS, while Tran- 
quada, et al.[17] showed the presence of nanoscale 
phase separation in the form of stripes, which ap- 
peared to be suspiciously coextensive with the pseu- 
dogap regime. After writing an extensive review[2], 
I put my ideas together into a self-consistent three- 
band slave boson calculation, reported at the first 
Stripes Conference[18]. Correlation effects keep the 
Fermi level close to the VHS from half filling to opti- 
mal doping. Near half filling, correlation effects drive 
the Cu-0 hopping to zero, and the remaining disper- 
sion due to J has a VHS at half filling, and gains 
additional stability by splitting the VHS via a flux 
phase[18,19]. Doping restores the hopping, simul- 
taneously introducing a strong electron-phonon cou- 
pling via modulation of the Cu-0 separation, leading 
to a maximal CDW instability at optimal doping. 
The two instabilities in turn drive phase separation. 
Attempts to model this phase separation led to good 
fits to the doping dependence of the photoemission 
dispersion. There was an important prediction: the 
VHS is found in photoemission to be below the Fermi 
level, because it is simultaneously above the Fermi 
level: the pseudogap consists of a splitting of the VHS 
into two features at (7r,0), but split in energy about 
the Fermi level. Photoemission could not reveal the 
upper VHS, but recent timneling studies[20,21] fully 
confirm this prediction, as well as demonstrating that 
the lower peak coincides with the photoemission VHS 
feature. 

At the mean-field level, the CDW has long-range 
order. However, when fluctuations are included in a 
mode-coupling scheme, there is only short-range or- 
der, with a real pseudogap opening up in the temper- 
ature range between the mean field transition tem- 
perature and a much lower transition to long-range 
order, driven by interlayer coupling[22]. If this inter- 
layer coupling is absent, the CDW resembles a quan- 
tum critical point (QCP), with correlation length di- 
verging as ^ ~ T~^/^[22]. It is not a conventional 
QCP, in that in the absence of phase sepatation it is 
not the terminus of a finite temperature phase tran- 
sition (i.e., there is no renormalized classical regime). 

2. SUPERCONDUCTIVITY and TUNNEL- 
ING 

The three-band model[18] revealed that the stri- 
ped phases are associated with two nesting instabili- 
ties, flux phase near half filling and CDW at optimal 
doping, leading to a large pseudogap near (tt, 0) and 
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simultaneously to nanoscale phase separation be- 
tween magnetic and hole-doped stripes. In the LSCO 
system, adding Nd or replacing the Sr with Ba leads 
to long range stripe phase order [17] while suppressing 
superconductivity - clearly demonstrating that three 
distinct phases must be involved. In the one-band 
model[12,15] this situation is simplified by replacing 
the striped phase by a uniform CDW phase. The key 
approximation is the Ansatz of reducing the strong 
correlation effects and phase separation by their most 
important effect: the pinning of the Fermi level to the 
VHS, over the full doping range from half filling to 
optimal doping. In the overdoped regime, the sim- 
plest approximation is assumed: the band structure 
stops changing with doping, and the Fermi level shifts 
in a rigid band fashion. 



Fig. 1. Phase diagram of pinned Balseiro-Falicov model. 
Circles = net tunneling gap, Aj. Inset: Tunneling spectra 
of a density-wave superconductor, using parameters of dotted 
line in main frame. Temperatures (from top to bottom) = 
130, 110, 90, 80, 70, 50, and 30K (dashed lines: T above the 
superconducting transition temperature). 

The pseudogap seen in tunneling measurements 
on Bi2Sr2CaCu208+5 (Bi-2212) [20,21] can be fit to 
this model. We have calculated the tunneling spec- 
trum[23]. Figure 1, and find a smooth evolution from 
pseudogap (CDW) phase to mixed CDW-supercon- 
ductor. For a pure CDW, the spectral function is of 
BCS form: 

Aik,i^) = 27r[ulS{Lu-Ek+)+vlS{iO-Ek-)]. (1) 
with ul^l-vl^ (1 + ek-/Ek)/2, Ek± = {ek+ ± 



Efe ± e/c+Q and Ek 



where the nesting vector Q = (tt, tt), and the gap Gk 
and dispersion are defined in Rcfs. [12,24]. The 
model involves three gap parameters, two (Go and 
Gi) associated with CDW order, and one A with su- 



perconductivity. Figure 1 shows the calculated phase 
diagram and the net low-T tunneling gap, defined 
as half the peak-peak separation. The inset shows 
that in the mixed CDW-supcrconducting state a sin- 
gle gap evolves in the calculated tunneling density 
of states p{E) (except for phonon structure). The 
ratio of the total gap to the CDW/superconducting 
onset temperature is nearly doping independent, 
2At/kBT* ~ 4.1. Since we use the model of Bal- 
seiro and Falicov (BF)[24] to describe the underly- 
ing CDW-superconductivity competition, we refer to 
this as the pinned BF (pBF) model. At present, it in- 
volves s-wave superconductivity, but we are working 
on an SO (6) generalization, including d-wave super- 
conductivity. 

For tunneling along the c-axis into a two-dimen- 
sional material, the tunneling density of states (dos) 
is an average of the in-plane quasiparticle dos[25]. In 
this case, there is a one-to-one correspondence be- 
tween features in the electronic band dispersion, as 
measured in photoemission, and peaks in the tun- 
neling dos. Fig. 2. The main tunneling peaks (A) 
coincide with the split electronic energy dispersion 
near (tt, 0) and (0, tt) of the Brillouin zone - and 
hence with the corresponding photoemission peaks, 
as found experimentally [21]. In the present BCS-like 
model, the slight discontinuity at the phonon energy 
produces a large peak in the tunneling spectrum (D). 
Note that at (tt, 0) the CDW and superconducting 
gaps combine to form a single feature (A) in both 
photoemission and tunneling, whereas the gaps split 
into two separate features, B associated with super- 
conductivity and C with the CDW, near (7r/2,7r/2). 



Fig. 2. Dispersion (a) and tunneling dos (b) near 
the Fermi level in the pBF model in the presence of com- 
bined CDW and superconducting order. Parameters are 
Go =7.45meV, Gi =14.3meV, A =13.1meV; other parame- 
ters as in Ref. [23]. Dotted lines = ghost bands. 
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The presence of only a single, combined gap at (tt, 0) 
is a consequence of an underlying SO (6) symmetry 
of the VHS, as discussed in the following section. In 
turn, this explains the smooth evolution of the pseu- 
dogap into the superconducting gap, Fig. 1 insert, 
which therefore need not be taken as evidence for 
precursor pairing in the pseudogap phase. 

A most exciting possibility is that by comparing 
the photoemission and tunneling data, one should 
be able to experimentally measure the pinning of the 
Fermi level to the VHS. Indeed, Renner, et al.[20], un- 
aware of this prediction[2], have noted that 'the pseu- 
dogap is centered at the Fermi level in both under- 
and overdoped samples. It is therefore unlikely that 
the pseudogap results from a band structure effect.' 



Fig. 3. Splitting of the Fermi level from the VHS, mea- 
sured as the normalized difference between the tunneling and 
photoemission pseudogaps, (Ap^ — A.tu)/ ^TU- Solid line 
= theory in the absence of pinning; open circles = derived from 
data of Refs. [21,26]. Inset: tunneling dos at several dopings. 

To quantify the extent of pinning, we test the 
null hypothesis. We compare the presently available 
data with the expected tunneling spectra for a pure 
d-wave superconductor in the absence of pinning, as 
doping is varied and the Fermi level passes through 
the VHS, Fig. 3. From the inset, it can be seen that 
the two peaks in tunneling arc symmetric when the 
Fermi level is exactly at the VHS (optimal doping), 
while the peak on the photoemission (inverse photoe- 
mission) side is stronger for underdoped (overdoped) 
samples, in accord with experiment. However, for 
strong enough doping away from optimal, the peak 
should split, which is not seen. Since the supercon- 
ducting gap shifts off of (tt, 0), we take the difference 
between the tunneling gap I^tu a-nd the photoemis- 
sion gap /S.pE as an experimental measure of the 
splitting. Fig. 3. While there is considerable error, 
the splitting in overdoped samples is close to what 



is expected, hut the splitting is absent in underdoped 
samples - strong evidence for Van Hove pinning. 

3. SO(6) 

A one-dimensional metal is susceptible to a vari- 
ety of instabilities, including singlet or triplet super- 
conductivity, CDW's and spin-density waves. These 
instabilities can be organized group- theoretically [27], 
either on the basis of a symmetry group, or in 
terms of a larger spectrum-generating algebra (SGA), 
which contains the hamiltonian as a group element, 
and hence can be used to generate the full energy 
spectrum. A similar analysis can be applied to the 
Van Hove scenario, in terms of an approximate 80(6) 
symmetry group, or S0(8) SGA[28]. 

This SO (6) group contains as subgroups Zhang's 
S0(5)[29], Yang and Zhang's SO(4)[30], and Wen 
and Lee's SO(3)[31] (SU(2)). It includes two 6- 
dimensional superspins, which form an 'isospin' dou- 
blet [28]: one combines Zhang's 80(5) superspin (an- 
tiferromagnetism plus d-wavc superconductivity) and 
the flux phase; the other involves s-wave supercon- 
ductivity and a CDW (as in the pBF model) with 
an exotic spin current phase. There is a most inter- 
esting evolution of these groups from one dimension 
to two. Table I. Lin, et al.[32] analyzed the group 
structure of a two- leg ladder. They found an SO (8) 
symmetry group, which involves the SO (6) group as 
a subgroup, plus operators which are antisymmet- 
ric for k ^ —k. These latter operators are essential 
in the one-dimensional case, where the Fermi surface 
consists of two points ±kpx, but are irrelevant for the 
VHS's, which are on the Brillouin zone boundary, and 
hence do not couple to these operators. Table I illus- 
trates the evolution of a single superspin (Lin, et al.'s 
d-Mott state) from 1-D to 2-D. In this Table, SDW = 
spin-density wave; sc = superconductivity; and 
Pi2 are symbols introduced by Lin, et al.[32] for two 
of their SO (8) operators, called band spin difference 
and relative band chirality; and a dash indicates that 
a corresponding operator is lacking. Note that the 1- 
D CDW connects fcp and -fcp, while the 2-D CDW 
discussed above connects (tt, 0) and (0, tt). Note fur- 
ther that the 80(6) structure persists down to lad- 
Table I. Superspin vs. dimensionality 
1-D ladder 2-D 
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ders two cells wide, and hence should remain valid in 
describing the striped phases. 

In the 1-D metals, the SGA aspect is more fun- 
damental than the symmetry aspect. The various 
instabilities are usually not degenerate in energy, as 
required by a symmetry group. Instead, they are gov- 
erned by the allowed interaction terms, g's - hence 
the name g-ology - and the object of research is to 
derive the allowed phase diagram as a function of the 
possible g-values. In this case, the SGA is useful in 
cataloging the allowed instabilities [27]. The situa- 
tion is similar for the VHS. The one-band model is 
not itself symmetric under S0(6), but shows consid- 
erable signs of the underlying S0(8) SGA. Thus, the 
form of the phase diagram in Fig. 1 is generic of any 
competition between a nesting operator and a pair- 
ing operator, while the pseudogap at (tt, 0) has the 
simple form 

' (2) 



where the A^'s are the gaps associated with the 
twelve components of both superspin vectors [33]. 
When t' = 0, this vector addition of the gaps holds 
for the full Fermi surface. Note that for a symme- 
try group, all the A,'s in Eq. 2 would have equal 
magnitudes. 

The corresponding 2-D G-ology phase diagram 
can be worked out [28,34], in analogy with the 1-D 
case. In the Hubbard limit, the only interaction is 
the U term, and the phase diagram has a natural 
evolution from SDW at half filling to d-wave super- 
conductivity in the doped materials. However, this 
simple picture cannot account for the striped phases, 
which compete with superconductivity (e.g., at 1/8 
doping, where there is long-ranged stripe order, su- 
perconductivity is suppressed) . By adding a phonon- 
mediated effective electron-electron coupling, a CDW 
phase can be stabilized in the doped material, and 
competition between CDW and SDW generates a 
striped phase. 

4. POSTSCRIPT: HOW WIDE ARE THE 
CHARGED STRIPES? 

There was considerable discussion at the Confer- 
ence as to whether the charged stripes were one cell 
wide or two cells wide (a.k.a.: site-centered or bond- 
centered). This may seem like a trivial issue, but it 
can have profound consequences: the nature of the 
striped phase at optimal doping, and ultimately the 
relevance of the t-J model to high-Tg superconduc- 
tivity. 



The issues may be conveniently addressed by 
considering the White-Scalapino[35] model of the 1/8 
doped state. The overall magnetic and charge or- 
ders are consistent with experiment [17], but our main 
concern here is the actual charge distribution. The 
charge periodicity is four Cu02 cells across, with 
two magnetic cells having an average hole density 
of ^ 0.07 holes per Cu, and two hole-doped cells, 
each containing ~ 0.18 hole. Since the magnetic and 
hole-doped stripes have equal width, this is the high- 
est doping at whieh a m,agnetic dom,ain wall m,odel 
makes sense. This is particularly true, since there is 
a strong tendency for the magnetic stripes to contain 
an even number of cells since in that case, there can 
be a spin gap, as in an even-legged ladder [36]. Hence, 
as the doping is increased above 1/8, the magnetic 
cells can no longer shrink in width, and there must 
be a phase transition in the nature of the stripes. [If 
the hole-doped stripes were only one cell wide, this 
anomalous behavior would only arise near a doping 
of 1/4.] 



Fig. 4. Charged stripes in the tJ model[37], interpreted as 

constant hole density domains of variable width. Open eircles 
= average hole density per row; dashed lines = guides to the 
eye; solid lines = phase separation model. 

White and Scalapino[37] have indeed reported 
evidence for such a phase transition; here we would 
like to present a reinterpretation of their data. Fig. 4a 
of Ref. [37]. Since the stripes are sensitive to bound- 
ary conditions, we prefer to look at the average hole 
density along each row (parallel to the stripes) of 
the simulation. Fig. 4. The data (open circles) fall 
very close to the form expected for a phase separation 
model (solid lines), with the same average densities 
as at 1/8 filling, but now the magnetic stripes retain 
their minimum width, while the hole-doped stripes 
get wider. Since 0.18 is close to optimal doping, the 
optimally doped materials are likely to be character- 
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ized by a set of widely separated magnetic ladders, 
with little residual interaction. The physics will be 
dominated by the physics of the hole-doped stripes 
at their special doping. 

This is bad news for the tJ model. It was specif- 
ically designed as a highly simplified model, which 
retained just enough physics to accurately describe 
the cuprates near the insulating phase at half filling. 
It is highly unlikely that the neglect of the oxygens 
and electron-phonon interactions will continue to be 
valid in the new hole-doped phase. 

We would like to thank A.M. Gabovich for useful 
discussions about tunneling, and NATO for enabling 
him to visit us. Publication 743 of the Barnctt Insti- 
tute. 
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